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ABSTRACT 

Wc use a large dark matter simulation of a ACDM model to investigate the clustering 
and environmental dependence of the number of substructures in a halo. Focusing 
on redshift z = 1, we find that the halo occupation distribution is sensitive at the 
tens of percent level to the surrounding density and to a lesser extent to asymmetry 
of the surrounding density distrbution. Wc compute the autocorrelation function of 
halos as a function of occupation, building on the finding of Wechsler et al. (2006) 
and Gao & White (2007) that halos (at fixed mass) with more substructure are more 
clustered. We compute the relative bias as a function of occupation number at fixed 
mass, finding a strong relationship. At fixed mass, halos in the top 5% of occupation 
can have an autocorrelation function ~ 1.5 — 2 times higher than the mean. We also 
compute the bias as a function of halo mass, for fixed halo occupation. We find that 
for group and cluster sized halos, when the number of subhalos is held fixed, there is a 
strong anticorrclation between bias and halo mass. Such a relationship represents an 
additional challenge to the halo model. 

Key words: Cosmology: observations 



1 INTRODUCTION 

The halo model of galaxy clustering (see Cooray and Sheth 
2002 for a review) takes as one of its usual assumptions that 
the clustering of a dark matter halo is dictated solely by its 
mass. Observationally when looking at a cluster or group it 
is easiest to count galaxies (see e.g., Hao et al. 2010, Yang 
et al. 2005, for recent cluster and group catalogues), which 
are usually associated with the population of dark matter 
subhalos. The role that the number of subhalos plays in 
the clustering of halos and also how that number is affected 
by environment has the potential to be linked directly with 
observations. If halo mass and number of subhalos affect 
clustering independently then this is a challenge to a main 
assumption of the halo model. In this paper we examine the 
role of the number of subhalos in a halo on clustering in two 
different ways, first by studing the environments (nearby 
overdensity) of halos and second by computing the correla- 
tion function. 

Gao, Springel and White (2005) showed that the clus- 
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tering strength of dark matter halos depends not only on 
mass but also on formation time. This result used a high 
resolution large volume simulation (the Millennium Simu- 
lation of Springel et al. 2005) and was the first of a set of 
results showing that halo mass is only the first order driver 
of halo clustering. Subsequent studies confirmed the result 
(e.g., Zhu et al. 2006, Harker et al. 2006) and showed that 
other dependencies exist apart from halo mass, such as con- 
centration (e.g., Wechsler et al. 2006, hereafter W06), and 
halo spin and shape (Bett et al. 2007) . W06 also found a de- 
pendence on number of subhalos, and Gao and White (2007) 
(hereafter GW07) on the substructure fraction, both at fixed 
mass in that halos with more substructure are more clus- 
tered. It is these relationships that we will be investigating 
in more detail in this paper. 

The clustering of halos can be probed in a different, 
related fashion by measuring properties of the local envi- 
ronment. Again in this case to first order the dependence 
of halo properties on environment was shown to be lim- 
ited to mass by Lemson & Kaufmann (1999). Advances in 
simulation size and resolution enabled this type of analy- 
sis to be extended to measure smaller effects and different 
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statistics. For example, the abundance of halo substructure 
in a high resolution simulation of side length 21.4 /i _1 Mpc 
was shown by Ishiyama et al. (2008) to be dependent on lo- 
cal overdensity. Wang et al. (2011) find that the large-scale 
tidal field significantly affects all halo properties they stud- 
ied (assembly time, spin, shape and substructure) at fixed 
mass. White et al. (2010), concentrating on cluster-sized ha- 
los show that their environment influences and causes phys- 
ical correlations in many observational probes of mass and 
richness (e.g., subhalo number, lensing and velocity disper- 
sion). 

The analysis in this paper is an extension of some of this 
prior work, which has already shown in many ways that the 
subtleties of halo clustering are not explained by the sim- 
plest forms of the halo model. Here we use the language of 
the Halo Occupation Distribution (HOD, Berlind & Wein- 
berg 2002) to examine the environmental dependence of sub- 
structure. In looking at clustering, we examine samples of 
halos from the point of view of their bivariate distribution 
(mass and number of subhalos), which it turns out can lead 
to some interesting behaviours which may suggest further 
improvements in the halo model. 

The structure of the paper is as follows. In Section 2 we 
briefly describe the N-body simulation we use to study the 
subhalo population of halos along with their large-scale envi- 
roment. We give our measures of environment and describe 
our method for computing statistics that depend on halo oc- 
cupation at fixed halo mass. The environmental dependence 
of the HOD is examined in Section 3 and a resolution test is 
carried out. In Section 4 we examine the clustering of halos, 
measuring their large-scale bias as a function of mass and 
subhalo occupation. In Section 5 we summarize our results 
and discuss their implications for both the theoretical halo 
model and observational probes of dark matter halos. 



2 SIMULATION 

We have used P-GADGET, an upgraded version of GAD- 
GETS (see Springel 2005 for details of an earlier code ver- 
sion) which is being developed for upcoming petascale su- 
percomputer facilities, to run a large dark matter simulation 
of a ACDM cosmology. The cosmological parameters used 
were: amplitude of mass fluctuations, as — 0.8, spectra in- 
dex, n B = 0.96, cosmological constant parameter Qa = 0.74, 
and mass density parameter Q. m = 0.26. In the present 
work we use the simulation output at z = 1, the same as 
that analyzed in our previous paper using that simulation 
(Khandai et al. 2011). The initial conditions were generated 
with the Eisenstein and Hu (1998) power spectrum at an 
initial redshift of z = 159. The basic simulation parameters 
are: box side length, 400 /i _1 Mpc, the number of particles, 
2448 3 = 1.5 x 10 10 , the mass of a particle, 3.1 x 10 8 /i _1 Mo 
and the gravitational softening length, 6.5/i~ 1 kpc. For ref- 
erence, our simulation volume is roughly half that of the 
Millennium Simulation (Springel et al. 2005) but our mass 
resolution is about a factor of three better. 

2.1 Halos and Subhalos 

Halos and subhalos are identified on the fly as the simulation 
runs, the halos using a standard Friends-of-Friends (FOF) 



groupfinder with linking length 0.2 times the mean inter- 
particle separation. The halo masses quoted are the sum of 
masses of all particles in the FOF group. 

We use the SUBFIND code (Springel et al. 2001) to con- 
struct a subhalo catalogue and to measure the mass for every 
subhalo. Groups of particles are retained as a subhalo when 
they have at least 20 bound particles, which corresponds to 
a minimum group mass of Af su b = 6.3 x 10 9 hT 1 Mq. We will 
investigate the dependence of our results on the minimum 
cutoff mass in the halo and subhalo masses. We also carry 
out a resolution test in Section 3.2 

The number of FOF halos in our simulation with mass 
greater than 10 10 /i _1 M© is 15.3 million. These FOF halos 
contain a total of 18.6 million subhalos (we do not differen- 
tiate between central and satellite subhalos). Several of the 
trends we will be looking at are relatively subtle and the 
statistical power coming from the large number of subhalos 
and halos is needed to make them readily detectable. 

We note that other authors have investigated the depen- 
dence of substructure on other halo properties (e.g., Jeeson- 
Daniel et al. 2011, Skibba & Maccio 2011) and clustering on 
substructure (e.g., GW07, Wang et al. 2011). Substructure 
can be defined in many different ways (such as the mass frac- 
tion in subhalos above a certain mass) and also be counted 
either within the FOF halo or within a certain radius (GW07 
try r2oo). We use the definition of substructure appropriate 
to the usual definition of the HOD, the number of subhalos 
above a certain mass within the FOF halo. We expect (as 
shown by GW07 and Bett et al. 2007) that the complex na- 
ture of the relationship between detailed halo properties and 
clustering will mean that different definitions of substructure 
can yield different results. 



2.2 Environment 

As a measure of environment we have a number of choices. 
The most standard is the overdensity within a sphere of fixed 
radius. We take this, using a fixed radius of 5 /i -1 Mpc co- 
moving to be our fiducial measure. We also test other radii, 
as well as using instead the overdensity in the sphere cen- 
tered on a halo that encloses a fixed mass 10 15 /i _1 M©. This 
Lagrangian definition of overdensity has been used before 
in the context of simulations by for example Colberg & Di 
Matteo (2008), who investigated the relationship between 
supermassive blackholes and their environments. 

We note that the overdensity within either a fixed ra- 
dius or Lagrangian volume is expected to correlate strongly 
with halo mass (see for example the extensive study carried 
out by Haas et al. 2011). Haas et al. 2011 show that it is 
possible to formulate measures of environment that are in- 
dependent of halo mass. In the present paper we have the 
more limited goal of investigating the dependence of num- 
ber of subhalos (rather than mass) on environment. We can 
sidestep the mass-environment correlation by either plotting 
a statistic we are interested in as a function of halo mass 
directly (the HOD) for different environments, or else by 
plotting a statistic (the correlation function and its relative 
bias) as a function of number of subhalos for fixed mass. Ei- 
ther way we have removed the dependence on mass for the 
purposes of our analysis. 
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Figure 1. The locations of the top 5% of halos by occupation at fixed halo mass (see Section 2.3) (red points) and bottom 5% (blue 
circles). We show only halos with a mass > 10 12 h~ 1 Mq in a slice of thickness 40 h~ 1 Mpc. 



2.3 Halo occupation at fixed halo mass: definition 

Once we have a halo catalogue and have enumerated the 
subhalos per halo, we would like to break up our full sam- 
ple into high and low occupation halos, where "occupation" 
is simply the number of subhalos (we do not differentiate 
between central subhalos or satellites). 

In order to avoid dependence of properties (such as clus- 
tering) on halo mass, we would like to somehow make the 
sample be effectively at fixed halo mass. One possible ap- 
proach involves applying an upper and lower mass threshold 
to the halo sample and then ordering the remaining halos 
by occupation. There are two disadvantages with this. First, 
applying mass cuts in this way will reduce the number of 
halos available for study, perhaps dramatically if the mass 
window allowed is narrow. Second, if the mass window is 
not very narrow one could legitimately worry that depen- 
dence on halo mass will still creep in, as there is obviously a 
strong dependence of halo occupation and halo mass. This 
means that for such a mass window, the top say 10% of ha- 
los by occupation could be significantly more massive than 
the bottom 10%. 

We avoid both of these problems by instead breaking 
our full halo sample into a large number of narrow bins 
in mass. In each one of the mass bins, we order the ha- 
los by occupation number per unit mass, choosing the top 
5% of halos, top 10 % and so on. When making a sample 



of halos chosen by occupation we then take the required 
fraction from each mass bin, so that we are left with a sam- 
ple that spans the entire mass range, but that was chosen 
by occupation at fixed halo mass. We have tried varying 
the width of the mass bins, Alog 10 (M), finding that below 
Alog 10 (M) = 0.5 our results are independent of its value 
(we use Alog 10 (M) = 0.2). 

We make first use of our set of halos ranked by occu- 
pation at fixed halo mass in Section 2.4 below when plot- 
ting their spatial distribution. The main use for this type 
of subsample will however be in Section 4 when we examine 
clustering. It will be useful to us to also look at halos ranked 
by mass for a fixed number of subhalos. In this case we will 
use the same technique, breaking the sample into a large 
number of narrow bins, but this time in subhalo number. 



2.4 Spatial distribution 

As an example of our categorization of halos by occupation 
we show in Figure [1] the spatial distribution of high and 
low occupation halos. Here high occupation halos are the 
top 5% by number of subhalos per unit mass at fixed mass, 
using the definition given above (Section 2.3), and the low 
occupation are the bottom 5 %. For reference, the mean 
number of subhalos in the high occupation sample is 17.6 
and the mean mass per halo is 3.8 x 10 12 hT 1 M©. The mean 
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number of subhalos for the bottom 5% sample is 4.1 and the 
mean mass per halo is 3.9 x 10 12 /i -1 Mq. 

Looking at Figure [JJ we can immediately see that there 
is a difference in the way the two subsamples are tracing 
out structures. Lower occupation halos appear to outnum- 
ber the high occupation halos in the low density regions that 
fill most of space. In the densest areas, there are many more 
high occupation halos. We will see in the following Sections 
3.1 and 3.3 that this is borne out quantitatively by consid- 
ering the halo occupation distribution in different density 
environments and also measuring the correlation function 
for the two samples we are plotting here. 



3 HALO OCCUPATION DISTRIBUTION 

Many observational measurements of galaxy clustering can 
be reproduced by theories that include two main ingredi- 
ents, the clustering of dark matter halos, and a model for 
the number of galaxies in a halo. The Halo Occupation Dis- 
tribution (HOD, see e.g. Berlind & Weinberg 2002, Zheng 
et al. 2002, Kravtsov et al. 2004) is a way of formulating the 
latter. In the case were galaxies are identified with dark mat- 
ter subhalos, the HOD can be measured from a simulation 
by simply counting the number of subhalos as a function of 
halo mass. In the present paper we will only concern our- 
selves with the mean number of galaxies in a halo, leaving 
the form of the probability distribution for further work. We 
note that Boylan-Kolchin et al. (2010) have found that sub- 
halo abundances are not well described by Poisson statistics 
at low mass, but rather are dominated by intrinsic scatter. 

The simplest assumption that can be made about the 
HOD is that the mean number of subhalos, iV su b depends 
only on halo mass. This assumption has been shown to hold 
relatively well by Berlind et al. (2003) in (hydrodynamic) 
simulations and is consistent with the e.g., Yang et al. (2005) 
observational measurements of galaxy groups. It is known 
however that substructure fraction does depend on environ- 
ments in simulations (e.g., Wang et al. 2011). We will see 
how this translates into changes in the HOD, first by looking 
at the local overdensity around halos and then some other 
measures of the environment. 

3.1 Density dependence 

We measure the density within a radius of each halo center 
of mass. We then rank the halos by density and plot the 
HOD for the halos as a function of halo mass for different 
density subsamples. Results are shown in panels (a) and (b) 
of Figure [2] where we show the HOD for the top 5% by 
density, bottom 5% by density and for the whole sample. In 
the top panels of that figure we show the fractional difference 
from the HOD for the whole sample for the two extreme 
density bins. To compute the error bars on the fractional 
difference, the volume was split into octants, and a jacknife 
estimator (Bradley 1982) was used to compute the error 
on the mean from the standard deviation of the jacknife 
subsamples. 

We can see from Figure [2{a) that the number of sub- 
halos in a halo does depend on the local density, with the 
halos located in the densest (5%) of environments having a 
peak difference in halo occupation of ~ 10% more subhalos 



than the set of all halos. This difference is most prominent 
for halos of masses 10 12 - lO 13 /!." 1 ^/© and becomes zero at 
lower and higher halo mass. The difference is even larger 
for halos in underdense regions, which have less subhalos 
than the set of all halos by up to ~ 40%, a result which 
again depends on halo mass. Looking at more extreme ends 
of the iVsub distribution, the peak shifts to the right (e.g., 
~ 10 13 ' 5 /i _1 M Q for the top 1% by occupation. 

The environmental dependence continues out to larger 
radius, as can be see in Figure [2jb) where we use density 
measured with 10 h~ 1 Mpc to rank halos. A third method to 
measure the density is that within a Lagrangian volume with 
mass 10 15 h~ x Mq. This is shown in panel (c), where the sign 
of the effect is the same, although the amplitude for the low 
density environments is very different (and therefore the ef- 
fect is detected at a lower level of significance). This is likely 
to do with the fact that the lowest density environments 
in this Lagrangian picture are being measured out to very 
large radii and therefore diluting the effect. For example, the 
mean radius that encloses 1O 15 /i -1 M0 for the "bottom 1%" 
line is 23 h~ 1 Mpc. 

The dependence of halo properties on environment has 
been investigated by many authors (e.g., Ishiyama et al. 
2008, Wang et al. 2011). Recently Jeeson-Daniel et al. (2011) 
have shown that environment does not correlate with sub- 
structure mass fraction on a halo by halo basis. This appears 
to be at odds with what we find here, but the definition of 
environment used by Jeeson-Daniel et al. (2011) is different 
to ours, being chosen so that it is not dependent on halo 
mass. Also we count subhalos by number and not by mass 
fraction. Wang et al. (2011) on the other hand have found a 
relationship between local tidal field of a halo and substruc- 
ture. Jeeson-Daniel et al. and Skibba et al. (2011) find that 
concentration (closely related to age) is more fundamental 
in setting a range of halo properties than mass. 

The destruction of substructure over time in halos can 
explain the anticorrelation between age or concentration and 
substructure (Gao et al. 2004). One might therefore expect 
there to be less substructure in halos which formed earlier, 
and have a higher concentration. Our results could be ex- 
plained therefore if halos in higher density regions had later 
formation times and lower concentrations, but this is not 
the case in general (as pointed out by W06 and GW07), so 
that the situation is more complex. In W06 and Wetzel et al. 
(2007) it was shown that low concentration halos and late 
forming halos do preferentially reside in high density envi- 
ronments, but only provided the halo masses are M > M». 
Zentner (2007) pointed out that is the general dependence 
that would be expected from an excursion set theory analy- 
sis (see also further work by Dalai et al. 2008). These papers 
argued that high-concentration and early-forming halos are 
found in high density environments when M < M, because 
these halos have their growth quenched by the tidal fields 
of nearby large halos. At z = 1, M, ~ 10 11 Mq so this 
effect is more relevant than at z = where almost all halos 
of interest have M > M,. 

Gioccoli et al. (2010) have plotted the HOD for halos 
identified at different redshifts, finding a systematic shift 
over all masses to more subhalos at higher redshifts (and 
higher concentration) , which is consistent with the interpre- 
tation of our finding above (see also Kravtsov et al. 2004, 
and Zentner et al. 2005). We do however only find a differ- 
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M halo/ h ~ lM M halo/ h ~ lM M halo/ h ~ lM 

Figure 2. The halo occupation distribution as a function of environment. The black lines in the bottom panels show the HOD for all 
halos. The red lines show the HOD for the top 5% of halos ranked by the local density and the other colours for other percentiles, as 
listed in the legend. The local density is given by averaging in a sphere of radius 5 h~ 1 Mpc around each halo (left), 10 /i _1 Mpc (middle), 
or the sphere that contains a mass of lO 15 /^ 1 Mq (right). The top panels show the fractional differences from the HOD for all halos. 
The mass limit for a subhalo in all cases was 6 X l(Ph~ 1 Mq. We show jackknife error bars on the points (see text). 




M halo/ h ~ lM M halo/ h ~ lM M halo/ h ~ lM 

Figure 3. The halo occupation distribution as a function of environment for different lower limits on the subhalo mass (panels from 
left to right). As in Figure [2] the black lines in the bottom panels show the HOD for all halos and the red and blue lines show the HOD 
for the top 5% of halos ranked by the local density and bottom 5%, respectively. The local density is given by averaging in a sphere of 
radius 5 h _1 Mpc around each halo. The top panels show the fractional differences from the HOD for all halos. We show jackknife error 
bars on the points (see text). 
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ence in HOD over a limited range of masses, so that other 
effects are involved as well. 

One thing which is clear is that the halo model of galaxy 
clustering which assumes a fixed HOD for all enviroments 
can be wrong at up to the 40% level in the least dense en- 
vironments. How this dependence of HOD on environment 
affects statistics that are measured, such as the correlation 
function of galaxies is best addressed by computing them 
directly, which we do in Section 4. 

Our minimum mass to be counted as a subhalo is 
6 x W 9 h~ 1 M e (20 particles). In Figure[3] we show the effect 
on the environmental dependence (defined as density within 
5 /i _1 Mpc) of the HOD on changes in this parameter. In 
panels (a)-(c) we are varying the subhalo mass threshold by 
factors of 3,4 and 16. The number of subhalos at fixed mass 
obviously decreases sharply as the subhalo mass threshold is 
raised. We however do not see any noticeable change in the 
difference between subhalo number in low and high density 
environments. This is interesting as one may have thought 
that small mass halos would be more susceptible to destruc- 
tion and so there should be more of an effect. This also does 
not appear to be a resolution effect (we return to this below) . 

In both Figure [2] and Figure [3] we can see that the en- 
vironmental dependence of the HOD appears to be much 
smaller, or go away completely for very high mass host ha- 
los (mhaio ~ 10 13 h~ 1 M(r ) ). This appears to be in accord with 
Zentner et al. (2005) and W06, who found that for low occu- 
pation number halos the environmental trends are those in- 
duced by formation time (because for small halos the abun- 
dance of subhalos reveals a lot about the mass accretion 
history of the host halo.) For larger halos this relationship 
is not as strong and so one might expect the environmental 
trend to be less evident also. 



p/<p> in 5 h 'Mpc, M sub >lO n M 

top 25% by p/<p> . . high res 

bottom 25% by p/<p> qq low res 




10" 10 1 
lo §io( M h a io/ M 0) 

Figure 4. A resolution test of the effect of local environment on 
the HOD. We show the fractional difference between the number 
of subhalos in all halos and in the top 25% picked by local density 
(red) and bottom 25% picked by local density (blue) (similar to 
top panels of Figure [2). Results for our fiducial simulation arc 
shown as smallcd filled points and results for a simulation with 
8 times worse mass resolution and smaller box size (see Section 
3.2) are shown as open circles. 



Given that this does not occur despite the large difference 
in mass resolution is some evidence for the robustness of our 
results. 



3.2 Resolution test 

As a test of the effect of resolution on our results, we have 
run a simulation with an identical cosmology but at worse 
mass resolution, and in a smaller volume. The simulation 
had 512 particles in a box of side length 167 /i -1 Mpc, so 
that the mass per particle is 8 times larger than in our fidu- 
cial simulation (the volume is also 14 times smaller). We 
again identified subhalos using SUBFIND, and compared the 
HOD in both simulations where the subhalo mass threshold 
was 40 particles in the low resolution run and 320 in the 
fiducial simulation (the same mass in each). 

The results are shown in Figure [4l where we show the 
fractional difference between mean number of subhalos for 
all halos and the top and bottom 25% selected by mass, 
(similar as in the top panels of Figures [2] and [3] except here 
we show only the top and bottom quartiles because of the 
small number of halos in the low resolution simulation). We 
can see the pattern again, with more subhalos in denser re- 
gions and a rapid dropoff at high masses in the number of 
subhalos in underdense regions. There are not enough halos 
with masses above ~ 5 x 10 12 M@ in the smaller, low 
resolution run to be able to compare results, but below that 
mass there does not appear to be any systematic difference 
between the N su b values and those for our standard simula- 
tion. If numerical effects were responsible for the destruction 
of subhalos in high concentration regions one would expect 
there to be significantly different amounts of substructure. 



3.3 Other environmental factors 

Although local overdensity is often used interchangeably 
with environment, other local measures of the environment 
of halos have been shown to affect halo properties. For ex- 
ample, Wang et al. (2011) find that at fixed halo mass, halo 
properties depend strongly on the local tidal field. The sub- 
structure mass fraction in their simulations is particularly 
affected, although the nature of the dependence is complex. 
They suggest that halos in higher density (and higher tidal 
field regions) have more accretion and so more substructure. 

We investigate three measures of local environment 
apart from the density and look at their effect on the HOD. 
Again we use a fixed radius of 5 /i -1 Mpc to define the local 
region around each halo. We look at the asymmetry of the 
local mass distribution (defined as the offset of the center 
of mass and the center of the halo, divided by 5 /i -1 Mpc), 
the mean infall velocity of particles, and the mean angular 
momentum. Of these, the asymmetry has some similarities 
with the tidal field studied by Wang et al. (2011), and the 
infall velocity will be strongly related to the overdensity. 
The results are shown in Figure [5l where we again plot the 
HOD for the top and bottom 5% of halos chosen according 
to each statistic. Again our fiducial cutoff in subhalo mass 
is 6 x 10 9 /i _1 M©. 

We compute the mean (mass-weighted) infall velocity 
by summing the mass times the component of the velocity 
of each particle in the radial direction (towards the center 
of the halo) and then dividing by the total mass of parti- 
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Figure 5. The halo occupation distribution as a function of environment for different measures of environment (all measured from 
matter in a 5 h~ 1 Mpc sphere surrounding the halo). The panels from left to right show respectively the asymmetry of the matter 
distribution (quantified by the difference between the center of mass of the matter with 5 h~ 1 Mpc and the center of the halo), the mean 
infall velocity of matter, and the total angular momentum of matter. As in Figurc[2]the black lines in the bottom panels show the HOD 
for all halos and the red and blue lines show the HOD for the top 5% of halos ranked by the environmental measure in question and for 
bottom 5%, respectively. The top panels show the fractional differences from the HOD for all halos. 



cles within the (5 ft _1 Mpc) radius. As expected, the high 
mean infall velocity selected halos (middle panel) show the 
same sign of effect in overabundance of subhalos as seen 
in highly overdense regions. The size of the positive effect is 
again close to a maximum 10% difference from all halos. The 
negative effect, in halos of low infall velocity is about a max- 
imum of 15% below the result for all halos, somewhat less 
than the equivalent result for density (left panel of Figure [2] 
The asymmetry results are more instructive, with the most 
asymmetric local regions having fewer subhalos (by ~ 2%) 
than all halos, and the least asymmetric ~ 2% more. That 
this is a separate effect to the density effect which can be 
understood by the fact that the most asymmetric regions ac- 
tually have higher densities (because they have more neigh- 
bouring halos). Looking at the relationship between density 
and asymmetry (not plotted) we find for the mean asym- 
metry, A ~ 0.15 + 0.25 log p 5 h -i M pc where p 5 h -i Mpc is the 
density within 5 /i _1 Mpc. We also find that the mean value 
of the density p 5 h -i Mp c! f° r the halos with the top 5% of 
asymmetry A values is the same as for the top 7% of ha- 
los ranked by density p 5 ;,-im P c- If increased local density 
is causing the trend, one should therefore expect (based on 
Figure 2) to find a 10% increase in the number of subhalos 
in the most asymmetric regions, rather than a 2% decrease. 



There appears to be little dependence of halo substruc- 
ture on angular momentum of the surrounding region (3rd 
panel), although the error bars are large. 



4 CLUSTERING 

Given that the environment on at least 5 /i _1 Mpc and 10 
ft _1 Mpc scales noticeably affects the number of subhalos in 
a halo at fixed mass, a natural progression is a study of the 
clustering properties of halos. As described in Section 2.3 
we have constructed samples of halos with varying numbers 
of subhalos and fixed halo mass, as well as samples of fixed 
mean halo occupation but varying mass. In this section we 
will examine the clustering of these various samples. Pre- 
vious work on substructure and clustering includes that of 
W06, who computed the mark correlation function of ha- 
los, with A^ub as the mark, and G07 who used substructure 
mass fraction as a defining variable for samples of halos.. 

4.1 Autocorrelation function 

In Figure [SJ we show the autocorrelation function £(r) of 
halos for three different mass thresholds, as a function of 
separation, r. In each case, we have selected samples of fixed 
mean halo mass but with occupation in the top 5% of halos, 
and bottom 5% of halos. The results for the middle panel 
(lower mass limit of 10 12 h~ x Mq) therefore correspond to 
the halos plotted in Figure [T] We can see that in each panel 
the autocorrelation function of the high occupation halos is 
significantly enhanced with respect to that of all halos (also 
shown). The low occupation halos, on the other hand have 
a slightly lower amplitude of clustering. 

We examine this relative bias in the top panels of Figure 
[5] where we plot the relative bias with respect to all halos, 
as a function of scale, e.g. for high occupation halos, 

&HOD/all(T) = \/£to P 5%(0/Call(r). (1) 
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Figure 6. The autocorrelation function £(r) of dark matter halos at fixed mass (Section 2.3) for different halo occupations. We show 
different lower mass limits on the halo mass in the panels from left to right. In each lower frame, £(r) for all halos above the mass limit is 
shown as a solid line and dashed and dotted lines are used to denote £(r) for the top and bottom 5% of halos ranked by occupation. The 
top frames show the relative bias &HOD/all( r ) obtained by dividing the curves for the top and bottom 5% of halos ranked by occupation 
by the curve for all halos. 



We can see in Figure [6] that the relative bias is close to 
fiat as a function of r for low occupation halos, and for high 
occupation it is flat on large scales, r > 5 h~ Mpc for all 
mass bins. The effect of high occupation is more pronounced 
for halos of larger mass, with the > 1O 13 /i _1 M0 halos having 
a relative bias &HOD/aii ~ 1.4 on large scales (an amplitude 
of £ ~ twice that of all halos, and the low mass cutoff halos 
m > 10 11 ft -1 Mq having bHOD/aii ~ 1-2 on these scales. 

Just as we expected from the plot of halo positions, 
there is a substantial clustering difference between halos of 
different occupation at fixed mass. How exactly this bias is 
related to halo occupation and mass can be examined by 
working with measurements of the large scale (flat part) of 
the bias directly, which we do below. 

4.2 Large-scale bias vs. mass and occupation 

Because a focus of this paper is on the link between clus- 
tering and halo occupation and mass we compute bias for a 
grid of values of these two parameters. To measure the large- 
scale bias b between halos and dark matter we use £(r) data 
points with r > 5 /i~ 1 Mpc. We describe our jacknife esti- 
mator for computing the error on b below. In Figure [7] we 
show results for & on a colour scale, as a function of halo 
mass and subhalo number along the two axes. The plane of 
the figure is not filled, but the relatively large scatter about 
the mean relation between N svi b and halo mass means that 
there is data for a substantial fraction of parameter space 
that lies off this relation. For example, at low mass and low 
iVaub there are 10 12 /i _1 Mq halos with only 1 subhalo of mass 
greater than our threshold (6 x 10 9 h~ 1 MQ), as well as some 
halos of similar mass but with iV au b a factor of 15 larger. 

The most obvious trend in the b values shown in the plot 
is the increase in b along the mean relation between iV su b and 



Mhaio- Halos of larger mass tend to cluster more and also 
on average tend to have more subhalos. It is also clear at 
low masses and low subhalo numbers that increasing halo 
mass but keeping iV su b fixed (moving along rows from left 
to right) yields an increasing b, as does increasing N su ^ but 
keeping Mhaio fixed (moving up columns). It is interesting 
that from this it appears that iV su b and Mhaio independently 
control the amplitude of clustering. This holds in the bot- 
tom left of the plot, but when one moves to the top right a 
trend along rows or up columns is harder to see. By using an 
averaging technique (Section 4.3 see below) we will see that 
it turns out that at high occupation, b is completely inde- 
pendent of Mhaio, (no trend along rows in Figure [7]) whereas 
at high masses, b is still correlated with iV au b (still a trend 
up columns, at least for the mass range shown). 

Before examining these trends in & as a function of halo 
occupation at fixed mass, we will explore the variation of 
b with halo mass and iV su b. The steepness of these trends 
will allow us to place our later results in context. That halo 
bias is strongly influenced by halo mass is a result that is 
at the heart of many analyses of large scale stucture (e.g., 
Kaiser 1984, Mo & White 1996, Seljak & Warren 2004) In 
Figure[S]we show 6 (in this case the bias of halos with respect 
to the dark matter distribution) plotted against mean halo 
mass in various mass bins, finding qualitatively the same 
steep relationship seen in e.g., Figure 8 of Seljak and Warren 
(2004). 

We also plot b against the mean number of subhalos 
in a halo, finding a steeper relationship (in that b increases 
faster for a given increase in logjV su b). Because jV su b in- 
creases more slowly than Mhaio as Mhaio is increased (see 
e.g., Figure 0) this behaviour is expected. We show & as a 
function of N su h for two different subhalo mass thresholds, 
finding that b increases faster for a higher threshold mass. 
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Figure 7. Large-scale bias of halos with respect to the dark 
matter as a function of halo mass and halo occupation. The au- 
tocorrelation function £(r) was computed for subsamplcs of halos 
defined by A^h a i os and iV su b. The large-scale (a fit to points with 
r > 5 fc _1 Mpc) bias b for these subsamples was computed with 
respect to the autocorrelation of dark matter particles (see Sec- 
tion 4.1) and is shown as a colour scale on the plot. Regions which 
are black are areas of parameter space (Mhaios; N au b) for which 
there are no halos. 

This again is expected as for higher mass thresholds the 
HOD curves stay longer on the shallower initial part of the 
A^sub — Mhalo relationship (see e.g., Figure 

4.3 Large scale bias at fixed halo mass and at 
fixed occupation 

We take samples of halos which are chosen to either have 
fixed halo mass or fixed number of subhalos (as defined in 
Section 2.3), and compute how the large scale bias ^HOD/aii 
changes as we vary the other parameter. Our results are 
shown in Figure [9] where on the x — axis we show the per- 
centiles of the distribution. For example in the top panel a 
point at x = +99% shows &HOD/aii for the top 1% of ha- 
los by occupation, at fixed halo mass. Likewise a point at 
x = —99% shows ^HOD/aii for the bottom 1% of halos by 
occupation, at x = —50% for the bottom 50% and so on. 
Points at x = show results for all halos in that halo mass 
bin, and the bias 6hod /all is measured relative to this (as in 
equation 1). This is different from b computed with respect 
to the dark matter distribution, and ensures that all curves 
pass through &HOD/aii = 1 at a; = 0. 

In the top panel of Figure [9l the values of &hod /ail are 
shown for variations in the -/V au b percentile, at fixed halo 




1 10 100 
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Figure 8. Large scale bias between halos and dark matter as 
a function of either halo mass (large points) or halo occupation 
(small points). For each point a lower limit was applied to the 
sample of all halos and the position of the point on the a;-axis 
corresponds to the mean value (of either Mhaio or ^sub) °f halos 
above the cut. For iV su b we show two curves for two different 
lower mass cuts on the subhalo mass. 

mass. We show results for two different halo masses, and 
two different cutoffs in the mass of a subhalos. The mean 
halo mass {Mhaio) shown in the figure caption for the curves 
was computed by averaging over the masses of halos above a 
mass threshold. As required, the {Mhaio) are approximately 
constant for the different bins of AT su b percentile, and in all 
cases within 10% of the (Mhaio) value given. We can see 
that &HOD/aii does change significantly with A^b percentile, 
even though the mean halo mass is the same for all points 
along the curve. This is what we expect from looking at 
Figures [T] and [6] We can also see that the steepest change in 
b occurs for the larger mass halos. This analysis is equivalent 
to adding together the trend in b that one gets by moving 
up different columns in Figure [7] 

If we now turn to the results showing how 6 varies as 
Mhaio is changed for fixed Asub (bottom panel of Figure [9]), 
we can see that the situation is somewhat different. In this 
plot, the {N su h) values shown are computed from the average 
number of subhalos above a threshold in Asub , for each of the 
bins of Mhaio percentile. Here we see that for a low number 
of subhalos, (N su b) = 1.2, there is a strong trend of b with 
halo mass, one which is strong for both values of subhalo 
mass cutoff. However when we increase the A^b threshold 
to give a higher value of (N au b) ~ 8, we find that there is no 
dependence of &HOD/aii on Mhaio percentile. 

This independence of halo clustering and halo mass is 
interesting enough that we revisit it in Figure [101 where we 
now turn to plotting 6 with respect to dark matter halo clus- 
tering. The scale of the effects can then be judged compared 
to the overall trend of clustering increasing with greater halo 
mass and with greater N BU b- 

In order to judge the significance of the comparison, it is 
useful to put error bars on the points. We have done this by 
first splitting the simulation into octants and using a jack- 
knife estimator. Using random catalogs we have computed 
£(r) for the halo subsamples and for the dark matter parti- 
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Figure 9. Large-scale bias &HOD/all as a function of halo occu- 
pation percentile for fixed halo mass (top panel) and halo mass 
percentile for fixed halo occupation (bottom panel). The defini- 
tion of the percentile scale is given in Section 4.3. The &HOD/all 
values shown are the relative bias between the halos and the 
subsample of halos in the percentile of either halo occupation 
of mass. All curves therefore pass through &HOD/all = 1 at x = 
(see Section 4.3). In the top (bottom) panel we show results for 
two values of mean halo mass (occupation). For each of these we 
show what happens when the subhalo mass limit is changed (from 
6 x 10 9 h- 1 M Q to 2 x 10 la h- 1 M Q ). 
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Figure 10. Large-scale bias (wrt. dark matter) for halos of fixed 
mean mass but varying N suh (top panel) and halos of fixed N suh 
but varying M^aio (bottom panel). In the top panel for reference 
we show b vs. iV su b for all halos as a solid black line. Each set of 
other lines in the top panel is for a fixed mean halo mass (given 
in the legend) and shows how b varies with jV su b in that case. In 
the bottom panel, b vs. M^aio is shown as a solid black line. Each 
of the other lines in the bottom panel is for a fixed mean JV su b 
(given in the legend) and shows how b varies with M^aio- I n all 
cases we apply the usual lower mass limit on the subhalo mass of 
6 X 10 9 h- 1 M P) . 



cles, for the simulation volume minus each octant in turn. 
The value of b was computed for each jackknife subsample, 
with the error bar coming from their standard deviation. 

In Figure [TTJ1 as in Figure [5] the top panel shows & as a 
function of A?s U b for fixed halo mass, and the bottom panel 
show 6 as a function of Mhaio for fixed iV su b . We show results 
for 4 different halo masses in the top panel, and the x— axis 
is iVaub instead of A^b percentile. We can see that even for 
-Mhaio = 4 x 10 12 /i -1 Mq there is a trend of increasing b 
with increasing iV su b. As in Figure [5] for the smallest AT su b 
we notice that there is some flattening and even a small 
increase in b. For larger N su h the curves for varying N su b 
at fixed halo mass appear to track the black line relatively 
well, which shows how b varies with 7V su b for all halos. 

In the bottom panel, the results for fixed N BU b show 



a different behaviour. We see that for a sample with A su b 
threshold of 4 subhalos (with mass > 6 x 10 9 h^ 1 M@ ) there 
is no dependence of b on halo mass. To show how different 
this behaviour is from that usually seen in the halo mass- 
bias relation (e.g., Seljak & Warren 2004) we can compare 
to the black line in the bottom panel of Figure [10] which 
shows how b varies as a function of halo mass for all halos. 
The relative change in b for all halos is substantial (a change 
of 35%) over the mass range that is probed by the blue line 
(where there is no change in b). 

The behaviour of this relationship is examined for larger 
halos in Figure 1111 We show the same bias versus mass plot 
as in the bottom panel of Figure 1101 but this time for halos 
with a lower limit on A su b varying from 10 to 100. We can 
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Figure 11. Large-scale bias (wrt. dark matter) for halos of fixed 
mean N an ^, but varying Af^alo ■ This plot is identical in nature to 
the bottom panel of Figure [TOl but shows results for larger halos. 

see that for fixed subhalo number, the lack of dependence 
of bias on halo mass seen in Figure \W\ gradually turns into 
a strong anticorrelation of bias on halo mass. This happens 
between a threshold of iV su b = 10 and A BU b = 20, halos 
which have masses of ~ 10 13 /i _1 M Q and above. 

We have seen earlier (Figure [2]) that the environmen- 
tal dependence of the HOD declines for halos larger than 
~ 1O 13 /i _1 M . From our clustering results ( Figure ITT)) one 
might perhaps have expected a negative environmental de- 
pendence which is not seen in Figure O The situation is 
undoubtedly complex. 



5 SUMMARY AND DISCUSSION 
5.1 Summary 

Using a large, high resolution dark matter simulation output 
at redshift z — 1, we have investigated the role that the 
environment of halos plays on the number of subhalos they 
host, and we examined the relationship between clustering 
and subhalo number. We find that 

(1) At fixed halo mass, the number of subhalos in a halo 
is affected by local density, with overdense regions having 
more substructures over the whole range of halo masses, 
and underdense regions having less. This effect can be as 
large as 40% for the most underdense 5% of regions. 

(2) Our finding (1) is not significantly affected by the 
mass limit applied to subhalos, or (in a resolution test) by 
simulation resolution. 

(3) At a much smaller level, the asymmetry of the local 
density (quantified by a centroid shift) affects the subhalo 
number, so that the most asymmetric halos have fewer sub- 
halos, at the percent level. 

(4) The clustering of halos at fixed mass is strongly af- 
fected by the number of subhalos. This is true over the entire 
mass range tested, from 10 11 - 10 13 /i -1 M Q and for differ- 
ent values of subhalo mass cutoff. As with prior examples of 
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the dependence on clustering of variables other than mass 
(e.g. age, Gao et al. 2005, concentration, W06), it is thus 
straightforward to generate samples of halos which have the 
same mass but widely different clustering properties. 

(5) The clustering of halos at fixed number of subhalos 
is only positively dependent on mass for small halos. As we 
increase the size of halos, we find that for halos with more 
than 4 subhalos (for a subhalo mass limit of 6x 10 9 /i _1 M o ), 
at fixed number of subhalos there is no dependence of clus- 
tering strength on halo mass, and then for larger halos, of 
group and cluster size, we find that this turns into a strong 
anticorrelation of clustering amplitude and halo mass. This 
case, showing first an independence of clustering and halo 
mass and then an anticorrelation of the two is perhaps the 
bluntest expression yet of the shortcomings of the standard 
halo model when dealing with the subtleties of galaxy clus- 
tering. 

5.2 Discussion 

The most widely used version of halo model of galaxy clus- 
tering has at its heart some very simple premises. It is 
already well known however from simulation studies that 
many parameters apart from halo mass are actually at work 
in the translation of halo clustering into galaxy clustering, 
including the formation time dependence of clustering (Gao 
et al. 2005), substructure dependence and concentration de- 
pendence (W06). Looking for weak spots in the predictions 
of the halo model can in principle make the search for bet- 
ter models easier. It can also highlight areas in which for 
the present non-linear simulations can be a necessary tool 
in the attempt to make precise predictions of galaxy cluster- 
ing. The areas we have investigated in this paper may have 
some observational consequences (as we discuss below), but 
they also serve to highlight some special cases for which the 
halo model's prime assumptions are in direct contradiction 
with what is seen. These include a sample of halos of the 
same mass which have radically different clustering proper- 
ties based on an internal property (subhalo number) as well 
as other samples of halos with very different masses but the 
same clustering amplitude. 

The environmental dependence of the HOD is a rela- 
tively subtle effect, as can be seen by the fact that some 
previous smaller simulations (e.g. Berlind et al. 2003) were 
consistent with no dependence. Also, the observed correla- 
tion function of galaxies can be well modelled by an envi- 
ronmentally independent HOD (Zehavi et al. 2004, Zheng 
et al. 2009). There are however already signs that some ob- 
served galaxy statistics are not well predicted by the halo 
model, including higher order clustering of SDSS galaxies in 
low density environments (Berrier et al. 2011). 

As future galaxy surveys move well past the million red- 
shift regime (e.g., Schlegel et al. 2011), we can ask however 
how well we need to know the HOD and any extra envi- 
ronmentally dependent terms in order to carry out precision 
cosmology with an HOD approach to galaxy clustering (e.g., 
Zheng et al. 2002). Further work is needed to determine the 
effect on for example the correlation function of including 
an environmental term in the HOD. We have also not inves- 
tigated in this paper the effect of environment on the distri- 
bution of the number of halos at fixed mass (e.g., Kravtsov 
et al. 2004) which may be strongly affected. 
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It is also not clear how relevant and how strong the ef- 
fects that we have seen here with dark matter simulations 
are on the galaxies that form within the subhalos. One can 
clearly enumerate many possible environmental effects that 
rely on baryonic physics (e.g., luminosities of backsplash 
galaxies, Pimbblet 2011, etc.). which will further compli- 
cate the effect of environment on the HOD, and which may 
even have the opposite sign. 

The obvious example where the number of subhalos is 
used to define a set of objects are optical cluster catalogs. 
The environmental dependence of halo occupation is likely 
responsible for a fraction of the scatter in the mass-richness 
relation noted in optical cluster finders (see e.g. Rozo et al. 
2011 for different sources of scatter). 

The cause of the HOD enviromental effect is a com- 
plex issue. It has been shown by many authors that a wide 
range of variables are affected by halo environment such as 
concentation (e.g. Wang et al. 2011). internal halo proper- 
ties such as substructure and shape are nearly all correlated 
(Jeeson-Daniel et al. 2011, Skibba & Maccio 2011) and the 
mass function of subhalos itself correlates with halo concen- 
tration, formation time etc (Gao et al. 2011). The relative 
importance of mergers and smooth accretion in building up 
the mass of halo (Wang 2011) is likely to play a role in 
the enviromental dependence of the HOD. Fakhouri and Ma 

(2010) have shown that mergers dominate halo growth in 
overdense regions and diffuse accretion in voids. If they sur- 
vive the merger and accretion process we therefore expect 
subhalos to be more numerous in overdense regions (for ha- 
los of a given mass), as we have found (see also Wetzel et al. 
2007). The destruction of subhalos over time by intrahalo 
merging and stripping will make this relationship even more 
difficult to decipher. 

Incorporating environmental dependences into the halo 
model, making it more complex, but able to deal with phe- 
nomena such as those that have been demonstrated here is 
one avenue which can be pursued. Recently, Gil-Marin et al. 

(2011) have presented some first steps in this direction. 
Another source of potential uncertainty in the predic- 
tions of the halo model is the definition of halo mass. More 
et al. (2011) have shown that the mass of a FOF halo in 
a simulation depends on resolution. More et al. state that 
the influence of substructures (which depend on redshift and 
cosmology) on the FOF halo boundary, will make it difficult 
to model this effect in general. In our work, the relation- 
ship between the number of substructures and environment 
is likely also influenced by the effect that substructures have 
on the mass definition. 

Direct computation of the dependence of the correlation 
function of halos for different occupations shows more of the 
complex relationship between halo properties and clustering. 
The bias of halos of the same mass can vary widely depend- 
ing on their occupation. For example (from Figure HT))) . the 
lowest 25% of halos by occupation at fixed mass can have 
a bias which is 25% lower than that for all halos. This is 
similar to the effect seen by GW07 based on substructure 
fraction within the FOF group (although not within r2oo) 
and appears to be a stronger trend than that based on other 
properties at fixed mass, such as formation redshift (GW05), 
concentration (W06) or spin (GW07). 

In this paper we have focused on the clustering of ha- 
los of galaxy mass, and also only looked at z = 1. Further 



work is needed to explore the relationship between clustering 
and halo occupation in cluster size halos and those at lower 
redshift. If the same relationships hold, then this could have 
interesting consequences for the clustering of galaxy clusters 
selected in optical surveys. One could make measurements of 
the dependence of galaxy cluster bias on mass (measured us- 
ing velocity dispersion or lensing mass) and richness, equiv- 
alent to halo occupation. Mapping out the bivariate distri- 
bution of bias values as in Figure [7] would give further clues 
to the nature of galaxy formation in groups and clusters 
and how much it is affected by the non-baryonic processes 
investigated here. 

Our final perhaps suprising finding is that one can eas- 
ily select samples of halos (by picking a fixed occupation) for 
which there is either no dependence of clustering on mass or 
even a strong anticorrelation between the two. This finding 
is one which could also be tested with observational data 
on both mass and occupation. It again points to the com- 
plexity of halo clustering and the difficulties of using galaxy 
clustering measurements for precision cosmology. 
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